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ABSTRACT 

Sersic parameters characterising the density profiles of remnants formed in collision-less disc 
galaxy mergers are obtained; no bulge is included in our simulations. For the luminous com- 
ponent we find that the Sersic index is n G (1.5,5.3) with (n) w 3 ± 1 and an effective 
radius of R e £ (1.6, 12.9) kpc with (R c ) « 5 ± 3kpc. The mean values of these quantities 
increases as the radial interval of fitting is reduced. A strong correlation of n with the central 
projected density Iq is found (n oc Iq 014: ) which is consistent with observations. No posi- 
tive linear correlation between the size (R c ) and structure (n) of our remnants is found; we 
do not advocate the existence of this. The photometric plane (PHP) of the luminous com- 
ponent (n oc i?" 05 /g 15 ) agrees well, within the uncertainties and the assumption of a con- 
stant mass-to-light ratio, with those observationally determined for ellipticals in the B-band 
(n oc R® 09 Iq 15 ) and for if-band remnants (n oc R® J§ ). We found that the surface 
defined by Sersic parameters {n, i? e , /io} in log-space is not a true plane, but a pseudo-plane 
with a small curvature at low values of n owed to intrinsic properties of the Sersic model. 
The dark haloes of the remnants have a 3-dimensional Sersic index of (n) « 4 ± 0.5 that 
are smaller than the ones obtained for dark haloes in ACDM cosmologies n « 6 ± 1. A 
tight dark Sersic "plane" (DSP) is also defined by the parameters of the remnants haloes with 
n oc ? , c ' 07 Po 10 ■ We conclude that collision-less merger remnants of pure disc galaxies have 
Sersic properties and correlations consistent with those of observed in early-type galaxies and 
local remnants. It seems that a "primordial" bulge in spirals is not a necessary condition to 
form bona fide ellipticals on grounds of the Sersic properties of remnants. 

Key words: galaxies: kinematics and dynamics - galaxies: formation - galaxies: fundamental 
parameters - galaxies: interactions - galaxies: elliptical - methods: TV -body simulations. 



1 INTRODUCTION 

Hierarchical galaxy formation theory (e.g. Cole et al. 2000, De Lu- 
cia et al. 2005, Bower et al. 2005) considers that early-type galax- 
ies have an accretion/merger origin, as was originally suggested by 
Toomre (1977). Observational (e.g. Schweizer 1998, Struck 2005, 
Rothberg & Joseph 2006, Kaviraj et al. 2006) and theoretical (e.g. 
Naab & Burkert 2003, Meza et al. 2005, Naab et al. 2005) evidence 
supports this picture although several topics remain unsolved (e.g. 
Peebles 2002, Tantalo & Chiosi 2004). 

Early-type galaxies show several correlations among their 
colours, luminosities, velocity dispersions, effective radii and sur- 
face brightness (e.g. Baum 1990, Faber & Jackson 1976, Ko- 
rmendy 1977, Djorgovski & Davis 1987, Dressier et al. 1987, 
Bemardi et al. 2003). These correlations provide constraints to any 
theory of formation and evolution of these galaxies. Furthermore, 
their properties are linked with the distribution of luminous and 
dark matter, that would be important when comparing with models 
of formation of elliptical galaxies. 

Observational studies [e.g. Caon, Capaccioli & D'Onofrio 
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1993 (CCD93), Graham & Colles 1997, Binggeli & Jerjen 1998, 
D'Onofrio 2001 (D01), Trujillo et al. 2004] have found that the 
surface brightness density profiles of early-type galaxies are bet- 
ter described by a Sersic (1968) i? 1,/n -profile than the classical 
de Vaucouleurs (1948) 7? 1//4 -profile. The index n is directly related 
with the curvature and "concentration" of the light profile (Trujillo, 
Graham & Caon 2001). 

Several observational relationships have been found between 
the index n and, for example, the total luminosity (L), effective 
radius (R c ) and central velocity dispersion (e.g. CCD93, Prug- 
niel & Simien 1997, Graham & Guzman 2003). Also, it has 
been found a linear relation among log n, log R c and /io (cen- 
tral brightness) termed the Photometric Plane (PHP) for early- 
type galaxies [e.g. Khosroshahi et al. 2000 (K00), Graham 2002], 
analogous to the Fundamental Plane (Djorgovski & Davis 1987, 
Dressier et al. 1987). Recently, Rothberg & Joseph (2004, RJ04) 
have found that nearby merger remnants have a peak in the n- 
distribution at n ~ 2 with most values in the range of 1 < n < 6, 
and in some cases it is found that n > 8. 

On other hand, theoretical studies of Sersic properties of 
merger remnants have appeared recently. For example, Gonzalez- 
Garci'a & Balcells (2005, GGB) and Naab & Trujillo (2005, NT) 
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find in collision-less simulations that bulge-less progenitors lead to 
ranges of n £ (2.4, 3.2) and (1.2, 3.1), respectively; when a sin- 
gle Sersic function is used to fit the entire remnant. For progenitors 
with a bulge component they obtain about the same range of Sersic 
index, n G (3, 8). Since "bona fide" ellipticals have values n > 4, 
they reach the conclusion that collision-less merger remnants of 
pure disc galaxies do not lead to concentrations, indicated by n, 
similar to those found in intermediate or giant elliptical galaxies 
(e.g. Graham et al. 1996). 

The above findings suggest that a primordial bulge in spirals 
is a necessary condition to form bona fide ellipticals in the hierar- 
chical merging scenario. However, we show below, collision-less 
mergers of pure discs can cover the range of observed values of the 
shape parameter n, and can reproduce adequately other observa- 
tional correlations. 

Sersic model in a de-projected form has been recently 
used to represent the dark matter distribution in ACDM haloes 
(Navarro et al. 2004, Merritt et al. 2005, Graham et al. 2005, 
Prada et al. 2005), in order to have a better estimation of the in- 
ner asymptotic logarithmic derivative. A mean value of a 3D Sersic 
index ~ 6, with a scatter of w 1, has been found in these works. So 
it is of interest to determine the three-dimensional Sersic parame- 
ters that characterise our remnants. 

In this work, we study the structural properties of remnants 
as provided by fitting a Sersic profile to their luminous and dark 
mass distribution. The paper has been organised as follows: in 
J3| we present a summary of the properties of our progenitors, 
some details of the simulations performed, as well as some basic 
characteristics of Sersic profile; both projected and deprojected. 
In S|3|we present distributions and correlations, in two and three- 
dimensions, found among the different Sersic parameters for our 
remnants, and compare them with observations. Sersic properties 
of the dark haloes of the remnants are determined, some correla- 
tions presented, and compared with those obtained in cosmological 
simulations. Some final comments are given in i|4]and a summary 
of our conclusions. 



2 SIMULATIONS AND SERSIC FUNCTIONS 
2.1 Galaxy models 

The galaxy models used in this work have been already described 
in Aceves & Velazquez (2005) and follow the method outlined by 
Shen, Mo & Shu (2002) to obtain the global properties of the discs, 
once the haloes properties are known. Our numerical galaxies do 
not include a bulge-like component. The dark haloes follow a mod- 
ified NFW (Navarro, Frenk & White 1997) model with an expo- 
nential cutoff. The discs have a typical exponential density profile, 
and satisfy the Tully-Fisher relation at redshift 2 = 1; roughly a 
look-back time of 8 Gyr in a ACDM cosmology with Hubble pa- 
rameter h — 0.7. Only discs satisfying the Efstathiou, Lake & Ne- 
groponte (1982) stability criterion were used. 

In this work, an additional simulation to those reported in Ta- 
ble 1 of Aceves & Velazquez (2005) has been done. This is a merger 
from the resulting remnants of A/01 and A/05, label MM. All 
simulations were carried out using a parallel version of GADGET- 
1.1 code, a tree base code (Springel, Yoshida & White 2001), and 
evolved for ~ 8 Gyr with conservation of energy better than 0.25 
percent. 



2.2 Density Profiles 

We fit only Sersic profiles to our merger remnants; no bulge-disc 
decomposition is attempted since progenitors lack any bulge com- 
ponent. The Sersic surface luminous-mass density profile is given 
by 

E(/?) = E e^ (fl/ ^ )1/ " , (1) 

where R is the projected spherical radius, R c is the effective radius, 
n the index of the profile, b = b(n) ~ 2n — 0.324 (Ciotti & Bertin 
1999) and Eo the central surface density. Index n is associated with 
the curvature and the concentration of the profile (Trujillo, Graham 
& Caon 2001); n = 1 corresponds to an exponential profile while 
the classical de Vaucouleurs (1948) profile is obtained for n = 4. 
The accumulated projected luminous mass, AZl(-R), is given 

by 

M^R) = £ E(i?)d(rf) = 2 ™T(y) x Rl , (2) 

where a = 2n, x = 6(i?/i? c ) 1//n , and ■y(a,x) is the incomplete 
gamma function. The total projected luminosity mass is given by 

M L = |j£r(2n) E f? 2 , (3) 

being T(a) the complete gamma function. A summary of Sersic 
projected profile properties is given by Graham & Driver (2005). 
When comparing our simulations with observations we assume a 
constant mass-to-light ratio, so that E cx /; where / refers to the 
surface brightness. 

The three-dimensional (3D) Sersic profile is 

p(r) =p e y ' °> , (4) 

where r is the spatial radius, d ~ 3n — 1/3 + 0.005/n 2 (Gra- 
ham et al. 2005) such that r c is the half-mass spatial radius. The 
total mass is determined from 

Aft = ^r(3n) p rl . (5) 

Sersic parameters for the luminous component were computed 
along 400 different random line-of-sights. To each projection a 
circularly averaged density profile E(i?) was determined, and a 
Sersic profile Q fitted by x 2 - m i mm i sa tion using the Levenberg- 
Marquardt method (Press et al. 1992) to obtain {n, _R e , Eo}. Sersic 
parameters for dark haloes are obtained by a similar procedure, but 
using equation @. 

2.2.1 Fitting Range 

The fitting set of parameters depend on the methodology used to 
obtain them. In particular, there have been indications that these 
parameters depend on both the covered range of surface brightness 
range (e.g. Capaccioli, Caon & D'Onofrio 1992) and the spatial 
radial interval for fitting (e.g. Kelson et al. 2000). 

Also, the determination of fitted parameters degrades when 
the inner parts of a galaxy are not well considered. For example, 
the index n tends more to be a representation of the outer slope 
of the profile than of the curvature of the luminosity distribution 
(Graham et al. 1996). The treatment and quality of data has also an 
effect on the fitted parameters. For example, CCD93 obtain higher 
values of n for NGC 4406, NGC 4552 and NGC 1399 (14.9, 13.9, 
16.8) in comparison with D01 (6.5, 7.2, 6.1). 

We have considered two radial intervals for our fits in order to 
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Table 1. Physical properties of remnants 



ID 


Rh 


fly 


K 


M lum 
10 lu 


M tot 


2T 

Jw] 


2 
1 




[kpc] 


[kpc] 


[km/s] 


[M ] 


[M Q ] 






M01 


66.9 


156.1 


213.0 


10.00 


16.60 


0.99 


0.32 


M02 


29.8 


71.4 


108.5 


0.60 


1.95 


0.99 


0.46 


M03 


24.6 


56.4 


99.0 


0.54 


1.29 


0.99 


0.53 


M04 


41.6 


96.2 


132.8 


1.55 


3.98 


0.99 


0.74 


M05 


22.2 


48.6 


100.3 


0.83 


1.15 


1.00 


0.93 


M06 


27.3 


63.4 


96.9 


0.81 


1.41 


0.99 


0.87 


M07 


24.0 


55.3 


105.8 


1.02 


1.45 


0.99 


0.51 


M08 


33.8 


80.5 


92.5 


0.37 


1.62 


0.98 


0.97 


M09 


28.7 


66.1 


103.5 


1.41 


1.66 


0.99 


0.98 


M10 


33.3 


74.9 


110.8 


1.66 


2.19 


0.99 


0.70 


Mil 


32.4 


76.6 


178.2 


4.47 


5.62 


1.00 


0.14 


M12 


32.1 


74.9 


147.0 


2.39 


3.72 


1.01 


0.18 


MM 


68.2 


163.1 


216.7 


10.72 


17.78 


1.02 


0.07 
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Table 2. Mean parameters using radial range Ii 



ID 


n 




R c 


-Mo 


—Mr 


RMS 








[kpc] 


[M /kpc 2 ] 


[M Q ] 




M01 


4.3 ± 0. 


.4 


9.2 ± 1.9 


28.5 ± 0.5 


27.6 


0.10 


M02 


2.1 ± 0. 


.1 


2.7 ± 0.2 


23.3 ± 0.3 


24.4 


0.10 


M 03 


1.9 ± 0. 


.1 


2.5 ± 0.0 


23. U ±0.1 


24.3 


0.12 


M04 


2.8 ± 0. 


.1 


3.9 ± 0.2 


25.1 ± 0.2 


25.5 


0.14 


M05 


3.9 ±0. 


.2 


1.7 ±0.1 


28.4 ±0.4 


24.8 


0.10 


M06 


2.5 ±0. 


.1 


4.2 ±0.3 


23.6 ±0.4 


24.7 


0.12 


M07 


3.1 ±0. 


.2 


2.6 ±0.3 


26.1 ±0.5 


25.0 


0.10 


M08 


2.6 ±0. 


.1 


2.1 ±0.2 


24.5 ±0.3 


23.9 


0.12 


M09 


2.7 ±0. 


.2 


8.1 ±0.8 


23.2 ±0.6 


25.3 


0.20 


M10 


3.2 ±0. 


.2 


6.4 ±0.7 


24.8 ±0.6 


25.5 


0.14 


Mil 


1.6 ±0. 


.1 


9.3 ± 1.4 


22.0 ±0.5 


26.6 


0.08 


M12 


3.2 ±0. 


.2 


4.3 ±0.6 


26.0 ±0.6 


25.9 


0.19 


MM 


2.4 ±0. 


.1 


9.1 ±0.7 


24.5 ±0.1 


27.6 


0.19 



asses their effect on the Sersic parameters. The first radial interval, 
Ji, is taken from our numerical resolution value £i = 100 pc to the 
outer radius 7795, which encloses 95 percent of the projected lumi- 
nous mass and is determined directly from the simulations; thus, 
Ii = [£1,7795]. The second one, I2, uses another inner point at 
£f = lO^i, 1 and outer point at 7770; this enclosing 70 percent of 
the luminous mass. For each line-of-sight used, two uniform ran- 
dom numbers £ £ [£i, £f] and t] £ [7770, 7795] are generated that in 
turn define I2 = In the Appendix we discuss some effects the 

radial range of a fit has on the parameters estimated using synthetic 
models. 



3 RESULTS 

In this section we present the results of the fittings done, both "lu- 
minous" and dark, to the merger remnants, as well as several rela- 
tionships among them based in observational studies. 

Table \\\ lists different global physical properties of our rem- 
nants obtained directly from the iV-body simulations. Column (2) 
is the total half-mass radius Rh, (3) the virial radius R v , (4) the 
virial velocity V v , (5) the total luminous mass Mi um and (6) the to- 
tal bounded mass Mtot , and column (7) is the virial ratio at the end 
of the simulation. The last column (8) provides the ratio of the to- 
tal mass of the secondary to the primary galaxy in the simulations. 
The merger labelled as MM corresponds to the simulation where 
the resulting remnants of M01 and A/05 were merged together in 
a parabolic encounter. 

Tables[2]and[3]summarise the mean values of the fitted Sersic 
parameters {n, R c ,fio} ipo = —2.5 log So), the total "magnitude" 
(M T = -2.5 log M L ) and the RMS of the fit, for the different pro- 
jections for both radial intervals I\ and I2; respectively. Here A/l 
is determined from the fitted values using equation Standard 
deviations are listed for the Sersic parameters. The values of Ml 
determined from the fits agree very well with Mi um . 



1 For reference, in a ACDM cosmology with h = 0.7 we have that 1" = 
464 pc at the distance of the Coma cluster (2 = 0.023), 977 pc at z = 0.05, 
and 4.5 kpc at 2 = 0.3 . 



Table 3. Mean parameters using random radial range I2 



ID 




n 




Re 
[kpc] 


-«) 
[M G /kpc 2 ] 


-M T 
[M ] 


RMS 


M01 


5 


.8 ±1.3 


13. 


.2 ±5.5 


30.9 ± 2.0 


27.8 


0.05 


M02 


2 


.1 ±0.3 


2. 


4 ±0.3 


23.7 ±0.8 


24.4 


0.05 


M03 


2 


.2 ±0.4 


2. 


3 ±0.2 


23.9 ± 0.8 


24.3 


0.05 


M04 


3 


.3 ± 1.0 


3. 


1±0.5 


26.5 ± 2.1 


25.4 


0.07 


M05 


3. 


.2 ± 1.1 


2. 


1±0.6 


26.6 ± 2.1 


24.8 


0.04 


M06 


2. 


.7 ±0.4 


3. 


7 ±0.3 


24.3 ± 0.9 


24.7 


0.06 


M07 


2 


.6 ±0.7 


3. 


0± 1.0 


24.8 ± 1.4 


25.0 


0.06 


M08 


2. 


.1 ±0.3 


2. 


3 ±0.5 


23.3 ± 0.8 


23.9 


0.04 


M09 


2 


.9 ±0.5 


6. 


9 ±0.9 


23.8 ± 1.2 


25.3 


0.10 


M10 


3 


.5 ±0.6 


5. 


8 ±0.7 


25.7 ± 1.4 


25.5 


0.06 


Mil 


1 


.6 ±0.3 


9. 


6 ±1.9 


22.0 ±0.8 


26.6 


0.06 


M12 


3 


.2 ± 1.4 


6. 


3 ±4.0 


25.5 ± 2.4 


26.0 


0.16 


MM 


3. 


.2 ±0.6 


9. 


2± 1.5 


26.3 ± 1.0 


27.6 


0.07 



3.1 Luminous Distributions 

3.1.1 Shape parameter 

FigureQ(?op) shows the frequency distribution of n for a set of ob- 
servational data in optical wave bands (DO 1 , La Barbera et al. 2005) 
and in the near-infrared (K) band [La Barbera et al 2005, Raviku- 
mar et al. 2005 (R05)]. A total of 169 galaxies in the optical and 
156 in the K band were used here. The frequency distribution of 
41 merger remnants observed in the if-band by Rothberg & Joseph 
(2004) are also indicated as a shaded histogram. The mean and 
standard deviations of these data sets are indicated, as well as their 
median. 

In Figure Q (bottom) we show the distribution of n for our 
merger remnants using the radial fitting intervals I\ and 7/2 • The fre- 
quency distribution for our iV-body remnants peak at a value n ~ 3 
in both cases; although using I2 it shows a somewhat broader dis- 
tribution. For Ji it is found that n £ (1.5, 5.3) and for I2 that n £ 
(1.4,9.5). These values are in good agreement with those found 
in intermediate mass ellipticals (e.g. Graham & Guzman 2003, de 
Jong et al. 2004, Trujillo, Burkert & Bell 2004, Ellis et al. 2005), 
some brightest cluster galaxies (e.g. Graham et al. 1996), dwarf el- 
lipticals (e.g. Binggeli & Jerjen 1998, Young & Currie 2001), and 
the local merger remnants of RJ04. 

Our results using I\, and no bulge, are consistent with the val- 
ues found by NT and GGB for their models with a bulge in the pro- 
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<R.> = 3 9±4.7 3.7±3.8 3.4±3 4 



P5 




2 - MERGS 



I, I, 
5.1±2.9 5.4±3.i 

4.1 3.8 



1 1 I , I 

1 2 4 6 8 10 

n 

Figure 1. (Top) Frequency distributions N(n) for Sersic index n of a sam- 
ple of early-type galaxies in the optical (OPT, solid line) and near-infrared 
bands (NIR, dashed line). (Bottom) Frequency distribution of our TV-body 
remnants using both fitting intervals I\ and I2 ■ The mean values and their 
standard deviations are indicated; below these, the median is indicated, it- 
band merger remnants (MERGS) of Rothberg & Joseph (2004) correspond 
to the shaded histogram in both panels. 

genitors. Furthermore, using interval I2 lead to some values n ~ 9. 
This does not seem to be due to the methodology in the computa- 
tion of the surface density profiles. NT construct artificial images 
analogous to the observational procedure while GGB fit ellipses 
to isodensity contours, both considering a wide range in the radial 
fitting range, and obtaining similar ranges for n. It is likely that 
differences in the way models of the progenitors are set up be prob- 
ably one of the reasons behind the differences with our results; see 



3.1.2 Effective Radius 

Figure [2] (top) shows the observed frequency distribution of effec- 
tive radii R c for the data considered in iQ.l.ll and that correspond- 
ing to our remnants (bottom). For the fitting radial range Ii we ob- 
tain R e G (1.6, 12.9) kpc, and for I2 we have R e £ (1.6, 34.5) kpc. 
The average value of the observational data is about 4 kpc and for 



1 , 1 . ~ 

0.5 1 5 10 50 

R. [kpc] 

Figure 2. Frequency distributions of R c for observed early-type galaxies 
(top) and TV-body remnants (bottom). Mean values and dispersion, along 
with medians, are indicated as in FigureQ The same observational data sets 
considered in FigurejTlare used here. 



our remnants is about 5 kpc. It can be noticed that R c shows a larger 
dispersion of values than the index n depending on the fitting inter- 
val. This was also noticed by Binggeli & Jerjen (1998). 

Our remnants have a lower bound of R c ~ 1.5 kpc, while the 
observational data considered here can reach smaller values R e > 
0.5 kpc. We are not able to reproduce the small values of R a mainly 
because in our sample of initial conditions no pairs of small progen- 
itors were included. On other hand, values of R c > 10 kpc can be 
reproduced by our more massive remnants (A/01 and MM); see 
Tables|2|and|i] 

A unique comparison with the distribution ofR a values found 
by NT and GGB is not possible, since their models can be scaled 
to arbitrary physical units; a thing that is not possible here due to 
the way our disc galaxy progenitors were built up. Nonetheless, 
if we use the range of dimensionless values found by NT (1 < 
R e < 1.7) for systems classified as pure "bulges", and use a length 
unit of 3.5 kpc (i.e., the radial scale-length of the Milky Way) to 
transform their results to physical units, we find that both results 
are consistent. Also, we obtain qualitatively the same behaviour as 



© 2004 RAS, MNRAS OOO.HTHol 



Sersic Properties of Remnants 5 





1 i i i i 1 i i | i i i i 




1 n 
1 u 


e 


o 




_ r z =0.54 




b 


RMS = 0.155 * A 

U ° O A 
O A . 


A^ 

:y * 
• «" 


4 


8 01*4* ,<#f 






- A ■■>& o "J; 


a* o^ 1 




... o>o . , 


qMERG 


2 




. N-body _ 






n kR°- 2z 




A 






r 2 =0.15 


1 




RMS=0.147 


0.8 


r z =0.15 




RMS = 0.187 





10 



0.5 



5 10 
R e [kpc] 



50 



Figure 3. Sersic index n versus effective radius for the observational data 
of D01, R05, and RJ04. Genuine E galaxies from D01 are represented by 
a triangle (A), while those classified as of an uncertain type by a star (*). 
Merger remnants of RJ04 are open circles (o). Lines correspond to least- 
square fits. The values obtained, along with the coefficient of determination 
and RMS in log n of the fit are indicated. The average results for our JV- 
body remnants, using the radial range I2, are shown with solid dots (•) and 
points correspond to different projections of those numerical remnants. 



the one shown in their Figure 1 8 where a sharp cut at the lower-end 
of the distribution, as well as an extended tail at larger values. 

Considering the observational values of R c and those in our 
TV-body remnants, we can establish with confidence that the sim- 
ulations can reproduce quite well the observed range of values. 
Even some large values of R c found in giant ellipticals (e.g., Gra- 
ham et al. 1996) are reproduced. 



3.2 Luminous Correlations 

Several works (e.g. CCD93, D01, R05) have found a series of cor- 
relations among Sersic parameters in early-type galaxies. We now 
turn to study some of these and compare them with the properties 
of our numerical remnants. Firstly, we consider two-dimensional 
correlations, and then turn to consider the so called Photometric 
Plane (PHP) [e.g. K00]. 



3.2.1 Two Dimensional Correlations 

In the work of Caon et al. (1993) it was stated that a linear positive 
correlation between n and R c exists for early- type galaxies; they 
find that n oc R®' 52 for early-type galaxies in Virgo. A similar 
conclusion was reached by D'Onofrio, Capaccioli & Caon (1994) 
analysing galaxies in Fornax. Combining the data of both works 
one finds n oc R®' 50 with a Pearson's linear correlation coefficient 
r = 0.72. The statement of CCD93 that structure (as indicated by 
n) of an elliptical depends on its size R c has been supported by the 
analysis of Trujillo, Graham & Caon (2001). 

Figure|3| shows index n against R c for some of the data con- 
sidered here, as well as the values obtained for our disc galaxy 



1 

0.8 
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Figure 4. Sersic index versus fiQ from observational data of D01 (A,*), 
RJ04 (o), R05 (x). and TV-body remnants mean values (•) and their pro- 
jected ones (points). Lines are the resulting scaling relations obtained from 
least-square fits. The coefficient of determination (r 2 ) and the RMS are also 
indicated. For clarity, results from ./V-body remnants have been displaced 
by a constant value along the /jg-axis. 



merger remnants. A linear least-square fit to the data of D01 leads 
to n oc Re' 37 , and for RJ04 mergers n oc i?c' 26 ; with linear corre- 
lation coefficients r = 0.73 and 0.39, respectively. A similar fit to 
our remnants yields n oc R® 22 with r = 0.39. 

However, the observational data plotted in Figure [3] shows a 
large scatter around the assumed linear correlation; a fact already 
noticed by other authors (e.g. Trujillo et al. 2001). These fluctua- 
tions are quantified by considering the coefficient of determination 
(r 2 ) that measures the proportion of the variance of one variable 
that is predictable from the other (e.g. Ryan 1997). As indicated 
in Figure |3 the coefficients of determination are rather small, and 
the RMS of the fits are large, so its is not clear that a true linear 
correlation exists between log n and log R c . 

We notice also that the TV-body remnant M05, being the 
smallest one, has the second largest n value in our simulations. 
These results lead us to state that there is no linear positive cor- 
relation between the "structure" and size of an elliptical. It seems 
that the values of n and R c are restricted by some physical mecha- 
nism to a finite region of the Sersic parameter space; an option also 
indicated by Trujillo et al. (2001). 

On other hand, a stronger observational correlation has been 
found between n and the central brightness /xq in ellipticals (e.g. 
K00, Graham & Guzman 2003) and it appears to extend to dwarf 
ellipticals (e.g. Binggeli & Jerjen 1998, R05). In Figure|4|we plot 
these quantities for the observational data of D01, R05, RI04, and 
for comparison our TV-body remnants. A linear fit log n-iio to these 
data leads to n oc /# 17 for D01, n oc 7^' 14 for both RJ04 and the 
ellipticals in R05, and n oc Jo' 14 for our merger remnants. All 
fits have r 2 > 0.9 and an RMS < 0.07, that lead us to conclude 
that log n-fto is a true linear correlation, at least for the range of n 
values considered. 

As shown, the numerical remnants presented here are able to 
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Table 4. Photometric Plane coefficients 


ID 




—b 




RMSri 


K00 


0.17 ±0.03 


0.069 ± 0.007 


1.2 ±0.1 


0.04 


D01 


0.09 ± 0.02 


0.058 ± 0.004 


1.3 ±0.1 


0.06 


K04 


0.21 ± 0.09 


0.074 ±0.013 


1.7 ±0.3 


0.13 


R05E 


0.15 ± 0.02 


0.066 ± 0.003 


1.1 ±0.0 


0.04 


R()5dE 


0.16 ± 0.04 


0.082 ± 0.004 


1.6 ±0.1 


0.07 


RJ04 


0.11 ± 0.04 


0.054 ± 0.004 


1.0 ±0.0 


0.06 


M(Ji) 


0.05 ± 0.00 


0.057 ±0.001 


-1.0 ±0.0 


0.05 


M(/ 2 ) 


0.05 ± 0.01 


0.057 ±0.001 


-1.0 ±0.0 


0.05 



reproduce very well the log n-po correlation and its tightness. We 
recall that we have assumed a constant mass-to-light ratio to con- 
vert Eo to To in order to compare with observations. Hence, it ap- 
pears that at least in the range of masses of our remnants (Table0, 
there is no need to assume a dependence of the mass-to-light ratio 
dependence with mass (or luminosity) to reproduce the observa- 
tions. 



3.2.2 The Photometric "Plane" 

Several authors [e.g. K00, Graham 2002, Khosroshahi et al. 2004 
(K04), La Barbera et al. 2005, R05] have found that Sersic param- 
eters {n, R c ,/j,o} of early-type galaxies define a plane in log-space 
of the form 

log n = a log R c + bp + c, (6) 

which is termed the "photometric plane" (PHP). Some authors 
instead of fio use the mean effective brightness {/i) G (e.g. Gra- 
ham 2002, La Barbera et al. 2005). The different 2D correlations 
of il3.2.1l can be considered then projections of the PHP. 

We have computed, by a linear-square fit procedure, the co- 
efficients of this PHP for our remnants under the assumption of a 
constant mass-to-light ratio. Since we find that using (/i) c leads to 
about twice the RMS in log n than using fio, we restrict ourselves 
to an expression of the form {6j. 

In Table|4]we list the values of the coefficients of J6j found in 
the works of K00, K04, RJ04, the elliptical and dwarf ellipticals of 
R05, and those we obtain from the data of D01. Also shown are the 
coefficients found for our merger remnants for both fitting radial 
intervals Ii and M(7~i) and M(/2), respectively. In Figure|5|we 
plot the PHP from these data, using for illustrative purposes the 
values from the data of D01 to define the abscissa axis. 

The remnants' coefficient b in ^6), associated with po, is rather 
consistent with the observed ones, aside of those found in dEs 
(R05). The coefficient a, associated with R B , is less well repro- 
duced. This is not surprising taking into account the large disper- 
sion in the n-R c relation (see Figure As several authors have 
pointed out (e.g. K00, K04, R05) a slight curvature towards small 
values of n is observed, a feature that tends to be reproduced here 
by the effect of merger M 11 that has n ps 1.5. The RMS is sim- 
ilar for both the observational data and our simulations. The best 
overall agreement is obtained with the data of D01. 

We consider that our a and b values are rather consistent with 
the whole set of values listed in Table|4| and that the numerical rem- 
nants are able to reproduce the PHP. In <)4]we argue that the PHP 
is not really a plane, but a "pseudo-plane" with a small curvature at 
low values of n due to the intrinsic properties of Sersic model. 
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Figure 5. Representation of the photometric plane (PHP) for the same ob- 
servational data of Figure^ and for our numerical remnants. Symbols are 
as in Figurel4l The a>axis is defined here with the values obtained from D01 
data. Scalings n oc R"Iq are indicated, along with the RMS of the linear 
fit. ./V-body remnants have been displaced a constant term. 



Table 5. Dark haloes 3D Sersic fits 



ID 


n 




log Po 


RMS 






[kpc] 


[M Q /kpc 3 ] 




M01 


3.3 


86.7 


9.08 


0.05 


M02 


3.5 


28.4 


9.81 


0.03 


M03 


4.5 


26.5 


10.90 


0.04 


M04 


4.0 


42.1 


10.17 


0.08 


M05 


4.3 


20.4 


10.93 


0.06 


M06 


4.4 


25.2 


10.88 


0.05 


M07 


3.9 


23.0 


10.32 


0.05 


M08 


3.8 


28.9 


9.99 


0.05 


M09 


4.5 


29.0 


10.92 


0.07 


M10 


4.0 


31.9 


10.24 


0.12 


A/11 


3.7 


41.3 


10.03 


0.06 


M12 


3.3 


39.5 


9.38 


0.05 


MM 


3.0 


82.6 


8.78 


0.04 



3.3 Dark Haloes 

Dark haloes in cosmological simulations are started to be- 
ing described by a 3D Sersic function (Merritt et al. 2005, 
Prada et al. 2005, Graham et al. 2005), of the form 

p(r) = po exp[-d n (r/r ) 1/n ] (7) 

with r being the spatial radius, r c a 3D "effective radius", and d n ~ 
3n - 1/3 + 0.005/n 2 (Graham et al. 2005). It has been found that 
provides a better fit to dark haloes than the typical NFW or M99 
(Moore et al. 1999) density profiles. Sersic indices of about 6, with 
a scatter of ~ 1, are found for the cosmological dark haloes. 

We have fitted 3D Sersic profiles to the dark haloes of our 
remnants. The radial range of the fit was from the convergence ra- 
dius r c (e.g. Power et al. 2003) to the dynamical virial radius of 
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Figure 6. (Top) Fitted 3D Sersic profiles to me dark haloes of our rem- 
nants. (Bottom) Logarithmic residuals of the fits. 

the remnant (Table0. However, instead of using the orbital period 
at the r2oo radius to determine r c as done in cosmological simula- 
tions, we used the orbital period at the virial radius. 

In Figure[6]we display the fitted 3D Sersic profiles, along with 
their residuals. In Table|5|the values of the 3D Sersic parameters are 
listed for each remnant, as well as the corresponding RMS. It can 
be seen that the 3D Sersic profile Q is a very good representation 
of the density distribution up to the virial radius of each remnant. 
This is in concordance with the behaviour of the 3D Sersic profile 
for characterising cosmological haloes. 

The haloes of the iV-body remnants have a mean 3D Sersic 
index (n) = 3.9 ± 0.5; the uncertainty being the standard devia- 
tion. A value that is lower than that found for cosmological haloes. 
However, it is consistent with the mean value of the dark haloes of 
the progenitors (n) = 3.7 ± 0.3; as expected from the preservation 
of the cuspyness of dark haloes in mergers (e.g. Boylan-Kolchin 
& Ma 2004, Aceves & Velazquez 2006, Kazantzidis, Zentner & 
Kravtsov 2006). 

Differences in results are expected since the outer radius of 
the fits are not the same; we use a dynamical virial radius while 
for cosmological haloes the fits are done up to r2oo or further out 
(e.g. Prada et al. 2005). It should be noticed that Sersic fits Q by 
Graham et al. (2005) have values for r c (see their Table 1) in some 
cases larger than their cosmological virial radius; these last ones 
listed in Table 1 of Diemand, Moore & Stadel (2004). For exam- 
ple, their haloes G02 and G03 have r c = 391.4 kpc and 405.6 kpc, 
respectively, while their virial radii are 337 kpc and 299 kpc; halo 
B09 shows even a more larger discrepancy. Unfortunately, they do 
not provide their numerical half-mass radii to make a direct com- 
parison with the values r c they obtained. Also, Merritt et al. (2005) 
and Prada et al. (2005) do not provide the fitted values r c and the 
numerical half-mass radii. This makes uncertain any comparison of 
our results with these works. 

Figure shows different relations among the 3D Sersic pa- 
rameters for the haloes of our remnants: n-po and po~r c , and in 
analogy to the PHP we have constructed a 3D dark Sersic plane 
(DSP). We find, assuming a log-linear correlation, that n oc po' 08 
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Figure 7. Relations between different physical quantities of the dark matter 
haloes obtained from a Sersic de-projected profile. 

and po oc r~ 3 ' 10 with coefficients of determination 0.96 and 0.70, 
respectively. This indicates that n-po can be considered with con- 
fidence a true log-linear positive correlation, as its was in the cor- 
responding 2D case, but we do not deem on the same level po-r e . 
The DSP found for the remnants is 

logn= (0.07±0.02) logr c - (0.04±0.00)p:o - (0.48±0.07) (8) 

where po = —2.5 log po- This turns out to be a very tight correla- 
tion, with a coefficient of determination of 0.98 and RMS of 0.007, 
for the range of haloes masses considered in our simulations (see 
TableEJ. 



4 FINAL REMARKS AND CONCLUSIONS 

The results found here, as well as those of NT and GGB, show that 
the merger scenario is capable of reproducing the Sersic properties 
of observed elliptical galaxies. This work shows, however, that the 
presence of a "primordial" bulge in the progenitors is not necessary 
to satisfy, for example, the observed values of the shape parameter 
n; as was suggested by NT and GGB. 

It is likely that the different results with NT and GGB have 
their origin on the initial properties of the progenitors; in particu- 
lar, their dark matter distribution. We have used cuspy (NFW-type) 
dark haloes in contrast to those used by NT (pseudo-isothermal) 
and GGB (Lowered Evans) that have a constant density core. It is 
probable that the higher concentration of dark matter used here, af- 
fected the distribution of luminous matter in a way to increase the 
index n which is correlated with the luminous concentration (Tru- 
jillo et al. 2001). Other initial conditions of the progenitors and of 
the encounters such as energy and angular momentum, both intrin- 
sic and orbital, may have played a role in the final concentration 
of luminous matter of the remnants, as indicated by the index n. 
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Figure 8. Form factor f2(n) in equation |5} contributing to the observed 
curvature of the photometric "plane" at low values of n; for an ensemble 
of equal total-luminous projected mass galaxies. An analogous f3(n) for a 
3D distribution is also plotted. 



A systematic study of the way different dynamical elements deter- 
mine the Sersic index is, however, beyond the scope of the present 
work. 

We have shown that haloes of remnants define a tight dark 
Sersic plane (DSP) analogous to that of the luminous matter and 
with less dispersion. No indication of curvature is present, at dif- 
ference to what is noted in the PHP at low values of n. We argue 
here that this curvature is real and is related to an intrinsic property 
of a Sersic profile. Consider the expression for the total luminous 
matter associated with a Sersic profile J3j. This can be written in 
log-space as 

log L T = log n - 0.4 /j, + 2 log R c + log f 2 (n) (9) 

with the "form factor" / 2 (n) = 2nT(2n)/b 2n . A given set of 
galaxies with equal Lt and different Sersic parameters would de- 
fine an exact log-plane, except for the presence of the log /2(n) 
term. In the 3D Sersic function Q me analogous form factor is 
/3 (n) = 47rr (3n) / d 3n . These non-constant terms introduce a sys- 
tematic change in a PHP-like expression. The importance of the 
form factor is larger for values n < 1 and smaller for n > 1; 
both /2(n) and fs(ri) are shown in Figure|8| Thus, the form factor 
of the Sersic model determines the curvature observed in the PHP. 
This explains why no curvature is found by La Barbera et al. (2005) 
whose galaxies show n > 2, but this can be seen in dwarf ellipti- 
cals with several values of n < 1 (K04). Also, the DSP does not 
show such curvature since n > 3 (see Figure^. Furthermore, the 
dispersion about these "planes" is determined by the luminosity or 
dark mass range of the galaxy sample. 

It remains to study the central phase-space densities of the 
remnants, to see if they are consistent with the estimates for ellip- 
ticals (e.g. Carlberg 1986), and to analyse their kinematical prop- 
erties with those observed in elliptical galaxies. We plan to study 
these topics in a future work. In summary, our main conclusions are 
as follows: 

(i) Collision-less mergers of pure disc galaxies yield values and 
distributions of Sersic parameters consistent with those observed 
for bona fide ellipticals. The existence of a bulge in merging spirals 



does not appear to be a necessary condition on grounds of Sersic 
properties of the remnants. 

(ii) The suggested positive log-linear correlation between the 
size (R c ) and structure (n) in ellipticals is not supported. However, 
the strong log n-fio linear correlation found in observational stud- 
ies is supported by our merger simulations. On other hand, the PHP 
is fairly well reproduced. For these results a constant mass-to-light 
ratio is assumed. 

(iii) The final dark haloes of remnants show values of n ~ 4 
lower than those found in cosmological simulations n ~ 6. The 
difference may be attributable to the non-equivalence outer radius, 
where the dynamical virial radius was used in our case to carry 
out the fitting by a Sersic profile. Haloes define a tight Dark Sersic 
Plane (DSP) in three dimensions, with no indication of curvature at 
the level of the smaller n obtained. 

(iv) The curvature observed in the PHP at low values of n is an 
intrinsic manifestation of the properties of Sersic model, due to the 
presence of a non-constant term dependent of n. 
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APPENDIX 

We briefly discuss here the effect of the radial interval on the fitting 
process of a Sersic profile to a mass distribution. To do this we 
generate exact Sersic profiles with R c = 1, Lt = 1, and n = 2,4, 8. 
Also, random fractional errors < {1, 10, 20}% are introduced. The 
radial fitting interval is chosen as follows. A random inner point £ 
is selected from the interval [0.03, 0.5]R C while the outer radius, 
r), is randomly generated from the interval [Rro, R95]; where 7?7o 
and i?95 correspond to the radii containing 70% and 95% percent 
of the projected mass. These two points define our radial fitting 
interval [£,rj\. To corroborate the importance the importance of the 
underlying mass distribution we use also a Hernquist (1990) mass 
model with an without errors. 



Sersic Distribution 

Table|6|lists the mean values of {n, R c ,fio} obtained from fitting 
1000 Monte Carlo experiments for three Sersic models with errors 
as indicated above. Each line lists, in order of the ascending error 
introduced to the theoretical Sersic profile, the parameters recov- 
ered from the fit inside the random interval [£,77]. The standard 
deviation for each quantity is provided. 

These results show that the determination of Sersic parame- 
ters is very robust, for errors < 10%, against the size of the fitting 
region. As the error in the ideal Sersic distribution increases the 
dispersion grows. This is more clearly appreciated for index n. In 
the limit of zero error, even for a random radial fitting interval, 
the model parameters are recovered exactly. For this case, we con- 
clude that the radial fitting range does not has an important effect 
on Sersic fitted parameters. 
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Table 6. Sersic fits. 



^true 




p 


— Mo 




5 


2 


2.000 ±0.015 


1.000 ±0.003 


0.956 ± 


0.030 


C 4 




2.014 ± 0.162 


1.002 ± 0.031 


0.981 ± 


0.319 






2.061 ± 0.418 


1.009 ±0.082 


1.067 ± 


0.819 


3 


4 


4.001 ± 0.037 


1.000 ±0.004 


4.940 ± 


0.080 


2 
1.2 




4.041 ± 0.408 


1.001 ±0.043 


5.026 ± 


0.876 




4.222 ± 1.749 


1.010 ±0.139 


5.407 ± 


3.719 


1 


8 


8.002 ± 0.087 


0.999 ± 0.007 


13.261 ± 


0.196 


ai 

0.8 




8.119 ± 1.042 


0.998 ± 0.071 


13.521 ± 


2.316 




8.663 ± 4.776 


1.006 ±0.251 


14.719 ± 10.407 


0.6 



Hernquist Distribution 

We now consider that the case where the underlying mass distribu- 
tion follows a Hernquist model. Here, 7?hi denotes its theoretical 
projected half-light (mass) radius. It is found that in the fitting in- 
terval [0.03, 2.79]i?hi , the underlying Hernquist's profile is fitted 
by a Sersic profile with index n — 2.6 and R c = 0.82 in agree- 
ment with NT. For a radial fitting interval of [0.03, 14.5] Rhi we 
find that n = 3.67 and R c = 1.10. This indicates that the process 
of fitting a Sersic profile is far more sensitive when the underlying 
mass distribution does not follows a Sersic one. 

The above was already noted by Boylan-Kolchin, Ma & 
Quataert (2005), where a systematic change in Sersic parameters 
was found when trying to fit a Hernquist profile. In Figure |9|(/e/0 
we reproduce this systematic effect, for the Sersic index n, R e , 
and the mean effective surface brightness (J c ). If a random error 
< 10% is introduced the trend is preserved but a large dispersion 
results; especially as the inner radius of radial interval of the fit is 
increased. 

Figure [Tol shows the distribution of fitted values with no er- 
ror (solid line) and with a random error < 10% (dashed line) for 
an underlying Hernquist profile where the radial interval was ob- 
tained from £ 6 [0.06, 0.50]i?hi and r\ € [7770, T795]. The mean and 
standard deviations of the distribution are indicated. For compari- 
son, the histogram of values corresponding to a Sersic model with 
n = 4 with a random error < 10% is also shown (dotted line); see 
Table|6| second line in the entry for n — 4. 

From the above results, it follows that when the underlying 
mass distribution is not of a Sersic type, the fitted values have a 
rather large dispersion even in the presence of no error. In partic- 
ular, higher values of the index n are obtained for different radial 
ranges of the fit. In order to have a confident estimate of n, and 
other parameters, one has to sample rather deep inside and outside 
the luminous (mass) distribution; from about 0.1 to 6 Rhi- 

In practise, for example, sampling very near the centre of a 
galaxy may pose problems due to resolution effects. This is a par- 
ticular problem for observations of galaxies at different redshifts, 
using the same angular resolution but representing different physi- 
cal scales, and can lead to uncertain Sersic parameters. 
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